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Three-flavor QCD simulation with the 0(a)-improved Wilson fermion action is made employing 
an exact fermion algorithm developed for odd number of quark flavors. For the plaquette gauge 
action, an unexpected first-order phase transition is found in the strong coupling regime (/3 < 5.0) 
at relatively heavy quark masses (mps/mv ~ 0.74-0.87). Strong metastability persists on a large 
lattice of size 12 3 x 32, which indicates that the transition has a bulk nature. The phase gap becomes 
smaller toward weaker couplings and vanishes at (3 — 5.0, which corresponds to a lattice spacing a ~ 
O.f fm. The phase transition is not found if the improved gauge actions are employed. Our results 
imply that realistic simulations of QCD with three flavors of dynamical Wilson-type fermions at 
lattice spacings in the range a — O.f-0.2 fm require use of improved gauge actions. Possible origins 
of the phase transition is discussed. 

PACS numbers: 11.15.Ha,12.38.Gc,05.70.Fh 



I. INTRODUCTION 

Realistic simulations of the strong interaction through 
lattice QCD require dynamical treatment of up, down 
and strange quarks incorporating their pair creation and 
annihilation effects in the vacuum. While simulations 
with the dynamical up and down quarks have now be- 
come routine (for recent studies, see Refs. 0, 0, 0, 0), 
adding a dynamical strange quark is still in the develop- 
ment stage. This is primarily because no exact algorithm 
to treat odd number of flavors was known until recently. 
In fact with the conventional Hybrid Monte Carlo (HMC) 
algorithm p|, the number of flavors is limited to even. 
The R algorithm [(| can be used for any number of fla- 
vors, as applied in recent three- flavor simulations with 
the staggered quark action 0, @ • However, the results 
are subject to systematic errors due to a finite step size 
in the molecular dynamics evolution, and this has to be 
controlled by taking the limit of zero step size which is 
quite computer time consuming. 

Recently, however, several exact algorithms for an ar- 
bitrary number of flavors have been proposed both for the 
Wilson-type 0, and the staggered-type 0] fermion 
actions. They have been shown to work with realistic 
lattice volumes without much increase of computational 
cost compared to HMC |Tol fTT| . Thus, a practical bar- 
rier to perform realistic three- flavor QCD simulations has 
been eliminated. 

The next step toward realistic simulations is to explore 
the parameter space of the three-flavor lattice QCD to 
ensure that the system is free from lattice artifacts over a 
range of lattice spacing a and quark mass m q adequate for 



phenomenology. In practice this means finding a region 
corresponding to a ~ 0.1-0.2 fm and m q = (1/3 ~ 1) x 
m s , with m s the physical strange quark mass. 

For the Wilson-type fermion action, there is a related 
important theoretical issue to settle. For even number 
of flavors, evidence has been accumulated over the years 
that vanishing of pion mass at a critical hopping parame- 
ter K c (fl) (/? = 6/g 2 with g the gauge coupling ) is due to 
spontaneous breaking of flavor-parity rotational symme- 
try El ■ Whether this understanding holds for odd 
number of flavors still needs confirmation. It is there- 
fore mandatory to cover the entire parameter space of 
the lattice action for three-flavor QCD. 

In this work we numerically explore the parameter 
space (/?, k) of three-flavor lattice QCD with the 0(a)- 
improved Wilson fermion action 14]. The three quark 
flavors are assumed to be degenerate in mass. For the 
gauge action, we test both the plaquette and renormal- 
ization group (RG) improved gauge actions. The 0(a)- 
improvement coefficient csw is fixed to the one-loop per- 
turbative value as fully non-perturbative values for the 
relevant gauge actions were not available until after the 
present work was well in progress 0, 0, 0] . 

For the standard plaquette gauge action, we unexpect- 
edly find [3 

a clear evidence of the existence of an or- 
dered phase for large k separated by a first-order phase 
transition from the disordered phase at smaller n. The 
transition persists for large volumes, and hence it is a 
bulk phase transition. Contrary to the disordered phase, 
the ordered phase does not exhibits the standard features 
of the confining phase. For example, the pseudo-scalar- 
to-vector meson mass ratio is similar in the disordered 
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and ordered phases at around 0.8, but the lattice spac- 
ing is unnaturally large in the ordered phase. The first- 
order phase transition is observed in the strong coupling 
regime. The gap of physical quantities across the transi- 
tion diminishes toward weaker couplings, and appears to 
vanish at ~ 5.0, which corresponds to a lattice spacing 
a ~ O.f fin. 

These results indicate that there is a parameter region 
which is not smoothly connected to continuum three- 
flavor QCD. In other words, the continuum three-flavor 
QCD can be approached only if one uses lattices much 
finer than a ~ 0.1 fm, if one employs the plaquette gauge 
action. 

We find that the metastability signals disappear if one 
employs improved gauge actions such as the renormaliza- 
tion group (RG) improved action or 0(a 2 )-improved 
Liischer-Weisz action |2(j. With these actions, the con- 
tinuum extrapolation should be possible from a conven- 
tional range of lattice spacings a ~ 0.1 0.2 fm. 

The finding of a first-order phase transition for the pla- 
quette gauge action is quite unexpected. A large number 
of simulations carried out in the past for the quenched 
and two-flavor cases were consistent with the expecta- 
tion that the disordered (i.e., confining ) phase extends 
over < k < k c for any value of [3. Recently, however, 
while making a study of twisted mass QCD, Farchioni et 
al. reported in a two-flavor simulation with unimproved 
Wilson quark action and plaquette gauge action |2l| that 
there exists a first-order phase transition at = 5.2. 
They suggested that this phase transition can be under- 
stood as the alternative to the parity-flavor broken phase 
which was pointed out by Sharpe and Singleton |22j. 

It is possible that their finding and ours have a common 
origin. Another possible explanation for our first-order 
transition is that it is related to the first-order phase 
transition encountered in pure SU{2>) gauge theory in the 
extended coupling space ((3, (3a), where a characterizes 
the strength of adjoint representation [23l IziL l25j. Fur- 
ther work is needed for clarification of the origin of the 
first order-phase transition both for the two and three 
flavor cases. 

The rest of this paper is organized as follows. In Sec- 
tion^we introduce the lattice actions and simulation al- 
gorithms we employed. Section Hni describes our study of 
phase structure of three-flavor lattice QCD with the pla- 
quette gauge action, where the presence of the metastable 
states is discussed in detail. The phase structure analysis 
for the case of improved gauge actions is discussed in II VI 
In Section we discuss possible origin of the first-order 
transition and related phenomena encountered in past 
studies. A conclusion is given in IVII where we also dis- 
cuss the possibility of realistic simulations of three-flavor 
QCD. The first report of the present work was briefly 
made in Ref. [Ta 



II. LATTICE ACTION AND ALGORITHM 



The partition function we study is defined by 
Z = JvU (det[D ud ]) 2 (det[D s ])e- s ^ u l 
Here S g (U) is the gluon action given by 



S g {U) 



(1) 



(2) 



where Wi x i and Wi X 2 ar e the plaquette and rectangular 
Wilson loops, respectively. The summations are taken 
over all possible plaquettes and rectangles on the lattice. 
The coefficients Cq and C\ are determined as Cq = 1 — 8cj , 
with ci = 0, -0.331, or -1/(12(P) 1 / 2 ) for the stan- 
dard Wilson action, the RG-improved action [l^, and 
the 0(a 2 )-improved (Liischer-Weisz (LW)) action [20| . 
respectively, and (P) is the plaquette average introduced 
for a mean field improvement. 

The fcrmionic determinant (det[Z? u( ;]) 2 represents the 
contribution of degenerate up (u) and down (d) quarks 
whereas the strange (s) quark effect is given by (det[Z? s ]). 
In this work we mainly consider the (9(a)-improved 
Wilson-Dirac operator D q = 1+M q +T q (q — u, d, s) with 
M q the usual hopping term including the hopping param- 



eter K q and T q the (9(a)-improvement SW term \l£ 
explicit form of T q is given by 

T q = "CgwKjdpyJ 7 ii V , 



The 



(3) 



with T^y the clover-leaf-type field strength on the lattice. 
For the coefficient csw we employ the value determined 
by tadpole improved one-loop perturbation theory as 



f 



csw 



<p)3/4 



1 + L sw ^ 



(4) 



where 4w = 0.0159 2 



27, 281 for the standard Wilson 



gauge action. For improved gauge actions it becomes 

4w = °- 008 ( RG ) or O- 013 (LW) H. The average pla- 
quette (P) is calculated in pure gauge theory with the 
same value of 0. 

We employ an exact HMC-type algorithm for three- 
flavor QCD developed in Ref. [Hj • The u and d quark de- 
terminant (det[D u d\) 2 is estimated by the usual pseudo- 
fcrmion integral: 



(det^]) 2 = folpfaexp [-\D^<f> ud \ 2 ] 



(5) 



To represent (det[D s ]) in a similar manner, we approxi- 
mate the inverse of D s by the non-Hcrmitian Chebyshev 
polynomial [3(], of order 2n : 



l/D s ^P 2n {D s ) = ^c^-iy 

n 

: : ]J(D s -Z* (k) )(D s -Z m ), (6) 



fe=l 



3 



with Zk = 1 — exp(j 2nk/(2n + 1)) [33] and a reordering 
index j(k). We then rewrite 



[det[D s ]) = det[D s P 2n (D s )} 

x Jv<t>\V<t> s vq> [-\T n {D s 



2 ], (7) 

with T n (D s ) ee J2toMD s - 1)'(= m=i(Ds - z m ). 
Introducing a fictitious momentum P conjugate to the 
link variable U, the effective Hamiltonian for the 2+1- 
flavor QCD reads 



1 . 
H = -P- 
2 



+ S g (U) + \D-i<t> ud \ 2 + \T n (D s )<f> s \ 2 . (8) 



We take account of the correction factor det[D s P2 n (D s )] 
by the noisy Metropolis test. After a trial configu- 
ration U' is accepted by the usual HMC Metropolis 
test, we make another Metropolis test with the accep- 
tance probability P cor r[^ — ¥ U'] = min [l,e~ dS ] with 



dS = lAiUT'M^xl 2 - \X\ 2 - Here A = [D s P 2n (D s )} 



1/2 

and x is the Gaussian noise vector with an unit variance. 
For other details of the algorithm, see Ref. [Icj . 



III. PHASE STRUCTURE FOR THE 
PLAQUETTE GAUGE ACTION 

A. Thermal cycle analysis 

We study the phase structure of three-flavor lattice 
QCD for the plaquette gauge action (ci = 0) assuming 
flavor degeneracy k ee k uc i = k s . Rapid thermal cycles 
are performed in the (/3, k) parameter space on 4 3 x 16 
and 8 3 x 16 lattices. For a fixed value of (3 we start 
a HMC simulation at n = and increase k in unit of 
0.001 at every 200 HMC trajectories. This process is 
continued until we reach the point at which we encounter 
a large violation of HMC energy conservation satisfying 
dH > 100. Then, we reverse the process and decrease 
k until it reaches the point sufficiently far away from 
the turning point. This procedure is repeated at a fixed 
interval over a range of (3. 

In this global scan, we do not perform the global 
Metropolis test of HMC or the noisy Metropolis test for 
the correction factor det[D s P2n(D s )] in Eq. JJJ. The or- 
der of the Chebyshev polynomial is fixed to 2n = 100 
and the molecular dynamics step size to dr = 1/40 em- 
ploying t = 1 for the length of unit trajectory. The 
stopping criterion of the BiCGStab solver is such that 
the residual defined by \Dx — b\/\b\ becomes smaller than 
10~ 14 (10 -8 ) for Hamiltonian (force) calculation, where 
D is the even-odd preconditioned 0(a)-improved Wilson- 
Dirac operator, x is the solution vector, and b is a source 
vector. The expectation values of observables are mea- 
sured during the last 100 trajectories after 100 thermal- 
ization trajectories at each k in the cycles. 

In Figure^we present the plaquette expectation value 
(P) during thermal cycles on a 4 3 x 16 lattice. The value 
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FIG. 1: Thermal cycles of averaged plaquette (P) on a 4 3 x 16 
lattice at /3 = 4.6-5.6 with 0.05 steps and 5.8, 6.0 from bottom 
to top. Each points are measured on 100 trajectories followed 
by 100 thermalization trajectories from its previous point. 
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FIG. 2: Thermal cycles of averaged plaquette (P) on a 8 3 x 16 
lattice at (3 = 4.6, 4.8-5.3 with 0.05 steps, and 5.4-6.0 with 
0.1 steps (from bottom to top) for the standard Wilson gauge 
action. Each points are measured on 100 trajectories followed 
by 100 thermalization trajectories from its previous point. 



of f3 increases from (3 = 4.6 to 5.6 from bottom to top in 
units of 0.05. We observe strong indication of metasta- 
bility at 4.8 < j3 < 5.1. The system appears to jump 
from a disordered phase at smaller k to an ordered phase 
at larger n. 

On small lattices such as 4 3 x 16, the gap and metasta- 
bility might be attributed to a "finite-temperature phase 
transition" due to a small spatial size, corresponding 
to the thermal first-order transition observed for three- 
flavor QCD with the (unimproved) Wilson fermion ac- 
tion [3^] . To examine this possibility, we show (P) during 
the thermal cycles on a spatially larger 8 3 x 16 lattice in 
Figure[21 While the range of (3 indicative of metastability 
is shifted and reduced, we still observe a clear hysteresis 
loop at f3 = 4.95 and 5.0. 

For the two f3 values we show the number of BiCGStab 
iterations iVcc required for the fermion matrix inversion 
during the calculation of the Hamiltonian JSJ with a given 
stopping criterion 10 - in Figure|3| We find a large gap 
in -/Vcg at the point of metastability signals. In the dis- 
ordered phase Nqg is relatively small, whereas in the 
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FIG. 3: Number of BiCGStab steps Ncg required for the 
fermion matrix inversion during the HMC trajectories at /? = 
4.95 and 5.0 for the standard Wilson gauge action. 



►~ (amps) in L phase 




* am v in L phase 




* (amps) 2 in H phase 





amy in H phase 




m 




m 111 
e s • 



7.30 7.35 7.40 7.45 

UK 

FIG. 5: Chiral extrapolation of (amps) 2 and amy at (3 = 4.9 
on a 12 3 x 32 lattice. In 'H' phase we observe larger plaquette 
values (ordered phase), and in 'L' phase smaller plaquette 
values are observed (disordered phase). 
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FIG. 4: A typical example of the two-state signal on a 12 3 x 32 
lattice at j3 = 4.88 and csw = 2.15. The plaquette history is 
shown. 500-1000 trajectories are devoted to the thermaliza- 
tion from ordered/disordered configurations. 



ordered phase iVcG is about a factor 2-3 larger. As we 
discuss in the next subsection the quark mass is heavy 
at the point of metastability. Indeed we find tops /my ~ 
0.83-0.87. Therefore the large value of iVcG; which im- 
plies a small condition number of the Wilson-Dirac oper- 
ator, in the ordered phase is not attributed to a physically 
small quark mass. 

We suspect the metastability and gap to continue to- 
ward strong couplings below j3 — 4.9. We cannot confirm 
it, however, since the thermal cycles encounter a large 
Hamiltonian difference dH > 100, and hence are turned 
back, before finding signals of the ordered phase. For 
example, at (3 — 4.8 this occurs at k = 0.137. At this 
point Ncg — 55.4, which is still a relatively small num- 
ber as seen in Figure [3J The same is true for j3 = 4.6 
and k = 0.134, for which Nqg = 35.6. Simply reducing 
the molecular dynamics step size or increasing the poly- 
nomial order does not resolve the occurrence of a large 
value of dH. Further studies are needed to understand 
the possible continuation of the ordered phase toward 
strong couplings. 



B. Exact simulation on a larger lattice 

To establish the nature of the hysteresis observed in 
Figure |2 we perform exact simulations starting from 
both ordered and disordered configuration with fixed val- 
ues of f3 and k on a 12 3 x 32 lattice. The ordered config- 
urations are made at a larger k value, and the disordered 
configurations are generated in the quenched limit k = 
at the same j3 value. 

Simulation parameters, including (3, cswj k j the num- 
ber of molecular dynamics steps Nmu , and the order of 
the polynomial 2n, are summarized in Table [fl We re- 
mark that a polynomial order of only 2n ~ 30-200 are 
necessary in order to achieve a ~ 90% acceptance rate 
(Peon) for the noisy Metropolis test. 

Figure 0] shows a representative result which demon- 
strates clear two-state signals persisting over 1000 tra- 
jectories. We confirm that the hysteresis seen in rapid 
thermal cycles are not an artifact of our inexact simu- 
lations in which the HMC and global Metropolis tests 
are skipped. Similar two-state signals are observed us- 
ing the i?-algorithm on a 8 3 x 16 lattice. Our observa- 
tion strongly suggests the existence of a first-order phase 
transition separating the ordered and disordered phases. 

In Table we list results for pseudo-scalar meson mass 
tops and vector meson mass my in lattice units. We have 
also calculated the Sommer scale tq using the condition 
r^dV (r) / dr\ r =r — 1.65 on the static potential. The lat- 
tice spacing estimated from the phenomenological value 
r = 0.49 fm is listed in Table D 

In Figure 13 we plot (aTOps) 2 and amy at j3 = 4.9 as 
functions of 1/k. If we extrapolate the data in the disor- 
dered phase (solid symbols) linearly as is conventionally 
made, we find k c = 0.13703(29) and my(n c ) = 0.67(11) 
which translates into a^ 1 = 0.87(14) GeV as a rough es- 
timate of lattice spacing. These are quite natural values 
comparable to those encountered in quenched and two- 
flavor simulations. However, before reaching this point, 
the three-flavor system makes a transition into an ordered 
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TABLE I: Lattice parameters used for the 12 3 x 32 simulations. (Pace) : HMC acceptance rate, Nmo : the number of MD 
step to proceed a unit trajectory. The molecular dynamics step size is given by dr = 1/Nmu- (Pcott) : Global Metropolis 
test acceptance rate for correction factor, 2n : the order of the polynomial. The symbols in the 'phase' column denote : L) 
smaller plaquette value (disordered phase), H) larger plaquette value (ordered phase), M) signals are much unstable and long 
autocorrelation over 500 trajectories is observed. 
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phase in which hadron masses are drastically reduced. 

If one looks at the lattice spacing determined from ro 
in Table |U we see that a" 1 « 0.8 GeV in the disordered 
phase, which is consistent with the spectrum estimate 
from TO p above, while a~ « 2 GeV in the ordered phase 
is much larger. Furthermore pion mass squared does not 
seem to decrease in the ordered phase. We then suspect 
that physical results cannot be obtained with simulations 
in the ordered phase. 



C. Phase diagram 

In Figure we plot the location of the observed phase 
transition nx((3) in the (j3 1 k) plane for various lattice 
sizes. Open squares for a 4 3 x 16 lattice and open cir- 
cles for a 8 3 x 16 lattice show the point of hysteresis 
observed with the inexact algorithm, while filled circles 
and open up triangles correspond to the point where two- 
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FIG. 6: Phase diagram in the (/3, «) plane. Data for the 
phase transition points Kx are from the rapid thermal cycles 
on 4 3 x 16 (open squares) and 8 3 x 16 (open circles) lattices. 
Filled circles and open up triangles are from the exact sim- 
ulation on 8 3 x 16 and 12 3 x 32 lattices, respectively. The 
critical hopping parameter k c (open down triangles) is esti- 
mated using an extrapolation from smaller k, values. 
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FIG. 7: Plaquette expectation values in the two phases on 
the first-order phase transition line. Data from the 8 3 x 16 
and 12 3 x 32 lattices are plotted as a function of (3. 



in combination with the plaquette gauge action. To guar- 
antee a smooth extrapolation to the continuum limit, we 
should carry out simulations at coupling weaker than the 
termination point of the first-order transition at (3 w 5.0. 
As shown in Tabled the lattice spacing estimated from 
r = 0.49 fm at (3 = 5.0 is l/a = 1.53(2) GeV (k = 
0.1320), 2.06(6) GeV (k = 0.1330) and 2.58(4) GeV (k 
= 0.1338). The largest n value (k = 0.1338) still corre- 
sponds to a heavy quark (mps/mv ~ 0.83). Taking the 
significant k dependence of l/a into account, the lattice 
spacing in the chiral limit would be even larger, possibly 
greater than 3 GeV. A 2 fm lattice would then require a 
30 3 x 60 volume or larger. Large-scale simulations start- 
ing at such fine lattices and large lattice volumes are too 
computer time consuming even with high-end supercom- 
puters available at present. 



state signals such as in Figure are observed with the 
exact algorithm. 

The location of the phase transition line significantly 
moves when we increase the lattice size from 4 3 x 16 to 
8 3 x 16. However, it stays at the same place when the 
lattice size is further increased to 12 3 x 32, which strongly 
suggests that the first-order transition line persists in the 
infinite volume limit at zero temperature. In fact, the gap 
in the value of (P) does not significantly change from 
8 3 x 16 to 12 3 x 32, as shown in Figured 

Figure [7| also shows that the gap in the plaquette ex- 
pectation value decreases toward larger (3 and vanishes 
around (3 = 5.0. Since we observe no sign of second-order 
transition at weaker couplings (Figure [2J, the transition 
appears to terminate at (3 ~ 5.0. 

Toward the strong coupling regime, we expect the 
phase transition line to continue below the last point of 
kx at (3 = 4.88 in Figure Indications are that the 
gap tends to become large toward lower (3 values (see 
Figure Q) , while the pseudo-scalar-to- vector meson mass 
ratio on the gap is almost independent of (3. 

The critical hopping parameter n c drawn by downward 
triangles in Figure[S|represents a rough estimate based on 
the number of iterations of the BiCGStab solver, -/Vcg> 
in the calculation of Hamiltonian (Eq. ©) on a 8 3 x 
16 lattice. Nqq is sampled during the thermal cycles 
and the estimates are obtained by linearly extrapolating 
1/Nq G as a function of 1/k. As mentioned earlier the 
thermal cycles do not extend far toward large values of k 
at the (3 values below 4.88. Therefore the estimate of k c 
may have a large uncertainty. We also emphasize whether 
there actually exists the critical k c where the pion mass 
vanishes in the ordered phase is an open issue. In fact 
the pion mass measured in the ordered phase at (3 = 4.9 
shown in Figure [S] is almost constant as a function of k. 

D. Practical implications 

Our findings expose a serious practical problem on sim- 
ulations using the 0(a)-improved Wilson fermion action 



IV. PHASE STRUCTURE FOR IMPROVED 
GAUGE ACTIONS 

If the first-order transition observed for the plaquette 
gauge action is a lattice artifact, one may expect that 
it can be eliminated by improving the gauge action since 
scaling toward the continuum limit is much improved and 
the lattice artifacts are expected to be suppressed for 
these actions. 

Here we test two types of improved gauge actions, both 
of which are defined with J2J : one is the RG improved 
action yjj and the other is the O(o 2 )-improved (Luscher- 
Weisz (LW)) action In the LW case, the tadpole fac- 
tor (P) in Ci is self-consistently determined at k = for 
each (3. The clover term is also determined by tadpole- 
improved one-loop perturbation theory (Eq. Q) with 
(P) at K = 0. Figure [S] shows the results of the thermal 
cycles on a 8 3 x 16 lattice. The simulation conditions such 
as the values for 2n and dr, and skipping of the HMC and 
global Metropolis tests, etc. are the same as those with 
the unimproved gauge action. In contrast to the case of 
the plaquette gauge action, we do not observe any rem- 
nant of hysteresis loop. This indicates that there is no 
first-order phase transition with these improved actions. 



V. POSSIBLE ORIGIN OF THE FIRST-ORDER 
TRANSITION 

A. Bulk phase transition in the 
fundamental-adjoint coupling plane 

In the pure SU(3) lattice gauge theory having both 
fundamental and adjoint couplings (3 and (3a-, & bulk 
first-or der p hase transition exists in the strong coupling 
regime |23l l24l |25| . The transition line starts at the 
purely adjoint point {(3,13a) = (0,6.5(3)) and extends 
toward larger (3 with decreasing (3a- It terminates at 
(4.00(7), 2.06(8)) and never crosses the purely fundamen- 
tal line (3a — 0, so that the pure gauge lattice theory with 
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FIG. 8: Thermal cycles of plaquette expectation value (P) 
on a 8 3 x 16 lattice, (a) with the RG improved gauge action at 
j3 — 1.50-2.25 in steps of 0.05, and (b) with tadpole improved 
Symanzik gauge action at /3 = 3.20-4.70 in steps of 0.10 (from 
bottom to top). 



B. Breakdown of parity-flavor symmetry 

For two-flavor lattice QCD with the Wilson-type quark 
action, Sharpe and Singleton [2^ carried out an analysis 
of parity-flavor symmetry using chiral Lagrangian tech- 
niques. They pointed out that, depending on the sign of 
an 0(a 2 ) term, there may be a line of first-order phase 
transition along which pion mass is non-zero and the chi- 
ral condensate flips sign, rather than a pair of second- 
order transition lines along which pion mass vanishes and 
parity-flavor symmetry breaks down. 

It is straightforward but complicated to extend this 
type of analysis to the three flavor case. Since there 
are three 0(a 2 ) terms allowed in the chiral Lagrangian, 
rather than a single term for the two flavor case, pre- 
dictions are less definite. Nonetheless one may similarly 
expect that first-order transitions may occur depending 
on the coupling of the three terms. 

Singular phenomena have been observed with dynam- 
ical Wilson-type fermion simulations in a variety of con- 
texts, and our finding is one more of the list of such 
phenomena. While it is not clear at present if the above 
analysis offer an understanding of these phenomena, we 
attempt to discuss them for orientation of future studies. 



Two flavor case 



the fundamental representation is smoothly connected to 
the continuum limit (3 = oo. Near the critical endpoint 
the correlation function of the ++ glueball channel di- 
verges, but the scaling of other observables is not much 
affected ||. 

This phase transition shares many properties with the 
one we found in three-flavor lattice QCD. The first-order 
transition separates ordered and disordered phases (large 
and small plaquette expectation values, respectively), 
and their lattice spacing measured through the string 
tension is largely different. Furthermore, the transition 
has a bulk nature, i.e. it remains in the infinite volume 
limit. 

One may suspect that dynamical Wilson fermions ef- 
fectively induce the adjoint gauge coupling, which gives 
rise to the bulk transition. This possibility was ex- 
plored for two-flavor unimproved Wilson fermion action 
in Ref. [34[ . They measured the strength of the induced 
adjoint coupling on the dynamical configurations and 
found that it is slightly negative, as opposed to the ex- 
pectation. The clover term might enhance the induced 
adjoint coupling as it has a 1 x 1 Wilson loop structure, 
but it is to be confirmed either analytically or numeri- 
cally. 



Farchioni et al. [2l( recently reported a first-order 
phase transition for the plaquette gauge action and unim- 
proved Wilson quark action at /3 = 5.2 and k ~ 0.1715. 
They suggested that this phase transition may be under- 
stood within the Sharpe-Singleton analysis. 

In an old work on the finite temperature phase tran- 
sition of two-flavor QCD with the Wilson fermion ac- 
tion, an unexpected strong first-order transition was 
found . It was suggested that this is a bulk phase 
transition since the Polyakov loop does not jump at the 
transition point below (3 ~ 5.0 while the plaquette ex- 
pectation value shows strong metastability. 

It is plausible that the two findings refer to the same 
bulk first-order phase transition. Numerically, Ref. |35| 
found a metastability at (3 — 5.22 and k = 0.17, in a 
close proximity of that in Ref. |2l|, and the plaquette 
values in the two phases reported by the two studies are 
in agreement. 

It is not clear if such the first-order phase transition 
persists when the Wilson quark action is improved by 
the addition of the clover term. For nqn-perturbatively 
0(a)-improved Wilson fermion action |36(, most physi- 
cal observables, such as hadron masses and matrix ele- 
ments, measured in the past simulations 0, l37| do not 
show singular behavior. On the other hand, the mass of 
++ glueball is surprisingly lower than in the quenched 
case pjq . perhaps hinting at the presence of a nearby 
singularity in the coupling constant space. Also the lat- 
tice artifact in the measurement of the light quark mass 
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FIG. 9: Thermal cycles of the plaquette expectation value 
(P) for the plaquette gauge action and unimproved Wilson 
fermion on 4 3 x 16 (top) and on 8 3 x 16 (bottom) lattices. 
The (3 values in the region 4.6-6.0 are scanned. Lines show 
different f3 values in steps of 0.1 (from bottom to top). 



through the axial- Ward- Takahashi identity is found to be 
rather large |3^ |. 

We also note that the strong first-order transition for 
the unimproved Wilson quark action disappears if the 
gauge action is improved [4(J |41| ■ 
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FIG. 10: Plaquette history of the Nj — 3 exact simulation 
with the plaquette gauge action and the unimproved Wil- 
son fermion action at f3 = 5.0. The hopping parameter k is 
changed from 0.1710 to 0.1718 as indicated in the plot. 



lattice (bottom). 

The absence of the gap on the 8 3 x 16 lattice does 
not necessarily mean that the chiral and continuum limit 
can be smoothly reached, since the thermal cycle above 
does not cover the small quark mass region. In order to 
cover this region we carry out a dedicated run at = 5.0. 
Figure ^| shows the history of plaquette over the HMC 
trajectories starting from k — 0.1710, where the thermal 
cycle ended, up to 0.1718. We observe large fluctuations 
of plaquette at k = 0.1716 and 0.1718, which my be hint- 
ing at the possible presence of a phase transition. 

Clearly further work is needed to reach a comprehen- 
sive understanding of the phase structure of lattice QCD 
with Wilson-type fermion action for two and three flavor 
cases. 



VI. CONCLUSIONS 



2. Three-flavor case 

The report by Farchioni et al. of a first-order phase 
transition for the unimproved Wilson quark action and 
the plaquette gauge action raises the possibility that a 
similar first-order transition may be present for the three 
flavor case. Indeed, in previous finite-temperature stud- 
ies, a lar ge la ttice artifact was found for this action com- 
bination : at the point of finite-temperature transi- 
tion, the light quark mass measured from the axial- Ward- 
Takahashi identity jumps for /3 < 5, contrary to the ex- 
pectation that the Ward- Takahashi identity holds at any 
physical phases with an identical value for the measured 
quark mass for the same bare parameters. 

We have attempted an initial thermal cycle study on 
4 3 x 16 and 8 3 x 16 lattices with unimproved Wilson quark 
action, and the results are shown in Fig.[5]for the (3 values 
in the range 4.6-6.0. We observe a signature of phase 

transition on the 4 3 x 16 lattice (top) at = 5.1 5.2, 

while metastabilities are not apparent on the 8 3 x 16 



We have reported the existence of an unexpected phase 
transition in three-flavor lattice QCD with the 0(a)- 
improved Wilson fermion action. It appears in the strong 
coupling regime f3 < 5.0 if one uses the standard Wilson 
plaquette gauge action, while there is no indication of 
such a phase transition for improved gauge actions. The 
phase transition persists for large lattice volumes, and is 
likely a bulk phase transition. 

Our findings pose a serious practical problem on simu- 
lations using the 0(a)-improved Wilson fermion action in 
combination with the plaquette gauge action. To avoid 
unphysical effects of the bulk transition, one has to carry 
out simulations at couplings weaker than its end point, 
but the lattice spacing is already smaller than a ~ 0.1 fm 
there, necessitating large lattice volumes and hence large 
computing resources. 

This circumstance motivates us to employ the RG- 
improved gauge action, for which we observe no strange 
phase structure, for large-scale three-flavor simulations. 
A non-perturbative determination of the improvement 
coefficient csw f° r a full O(a) improvement has been 
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made using the Schrodinger functional method 0, IT5 
ITU , and preliminary results on the light hadron spectrum 
have already been presented in Ref. |42l ]. 

Finally, the recent report that the two-flavor system 
with the unimproved Wilson action also has a first-order 
transition, and that it may be understood with the con- 
text of the Sharpe-Singleton analysis on realizations of 
the parity-flavor broken phase raises an interesting prob- 
lem that we need to clarify for phenomenological appli- 
cations of full QCD simulations with Wilson-type quark 
actions. 
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